Temperature Anisotropy in a Shocked Plasma: Mirror-Mode 
Instabilities in the Heliosheath 
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ABSTRACT 

We show that temperature anisotropics induced at a shock can account for 
interplanetary and planetary bow shock observations. Shocked plasma with en- 
hanced plasma beta is preferentially unstable to the mirror mode instability 
downstream of a quasi-perpendicular shock and to the firehose instability down- 
stream of a quasi-parallel shock, consistent with magnetic fluctuations observed 
downstream of a large variety of shocks. Our theoretical analysis of the solar 
wind termination shock suggests that the magnetic holes observed by Voyager 1 
in the heliosheath are produced by the mirror mode instability. The results are 
also of astrophysical interest, providing an energy source for plasma heating. 

Subject headings: instabilities — shock waves — solar wind 



1. Introduction 

Planetary bow shocks and interplanetary shocks serve as a unique laboratory for the 
study of shock waves in collisionless plasmas. Observations of these shocks usually show that 
ion distributions are anisotropic with respect to the background magnetic field downstream 
of the shocks. Mirror mode waves associated with this anisotropy are observed downstream 
of quasi-perpendicular shocks (defined by the angle between the shock no rmal and the up 
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dominates (e.g. JAnderson et al.lll994j ; ICzaykowska et al.ll200ll ). Magnetic fluctuations down- 



stream of quasi-parallel shocks {9bu < 45°) have not been identified in detail, but hybrid 
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simulations show that the firehose instability can occur downstream of these shocks for cer- 



tain ranges of upstream A 
at ~ 0.5 fcan fc Swift! 



fven Mach number (Ma) and plasm a beta, for instance, Ma > 3 



1983 



firm this point (e.g., iGreenstadt &: Fredricks 



1979; 



Krauss-Varban fc Omidil 1991); observations seem to con- 



Bavassano Cattaneo et al.ll2000l ). Quasi- 



perpendicular shocks are characterized by a sharp increase in the magnetic field strength, but 
quasi-parallel shocks are often more turbulent, with shock ramps containing large- amplitude 
waves which spread upstream and downstream. 



The recent crossing of the termination shock (TS) by Voyag er 1 (VI) (jBurlaga et al. 



20051 ; iDecker et al.ll2005l ; iGurnett fc Kurthll2005l ; IStone et al.ll2005l ) provides an opportunity 
to study shocks and shock-induced waves in the heliosheath. The sharp increase in the mag- 
netic field strength across the TS an d the downstream m agnetic field configuration suggest 
that the TS is quasi-perpendicular (jBurlaga et al.ll2005l ). As in planetary magnetosheaths 
downstream of a quasi-perpendicular bow s hock, the heliosheath s hows compressive magnetic 
fluctuations in the form of magnetic holes ( Burlaga et al. 2006al lbl). 



We propose a theoretical explanation for the temperature anisotropies and associated 
instabilities induced at a shock. Based on the theoretical analysis, we show that the magnetic 
holes observed in the heliosheath could be mirror mode fluctuations. The present results 
also provide a prototype for understanding shocks in various astrophysical contexts, such as 
gamma-ray bursts, supernova explosions and active galactic nuclei. 



2. Theory 



An anisotropic ion distribution requires the use of a pressure tensor in the magnetohy- 
drodynamic (MHD) equations. Observations show that the shock structure and dynamics 
depend on the shock geometry and upstream Mai and /3i, so we examine the temperature 
anisotropy A 2 = T± 2 /Tp as a function of these parameters, where T± and Tji are the plasma 
temperature perpendicular and parallel to the magnetic field. The subscripts 1 and 2 indi- 
cate physical parameters upstream and downstream of a shock. Assum ing a bi-Maxwellian 
plasma, we write the jump conditions across a shock as (IHudsonl 119701 ) 
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where /io is the permeability of vacuum, p is the plasma density, and v and B are the 
plasma velocity and magnetic field with subscripts t and n denoting the tangential and 
normal components with respect to the shock surface. The velocity is measured in the shock 
frame. The square brackets indicate the difference between the preshock (1) and postshock 
(2) states. The perpendicular and parallel plasma pressures are defined as P± = pk^T^jm v 
and P\\ = pk^T\\jm v ^ where &b and m p represent the Boltzmann constant and proton mass, 
respectively. For simplicity, we assume that the bulk velocity is parallel to the shock normal. 
The components of the magnetic field upstream of the shock are given by B nl = Bi cos 6 Bn 
and B t \ = BismOsn- Differe nt forms of the solutions to these equ ations have been derived 
based on various assumptions (IChao et al.lll995uErkaev et al.ll2000l ). The focus of the present 
analysis is to investigate under what conditions the shocked plasma is unstable to the thermal 
anisotropy instabilities. 

The thermal anisotropy serves a free energy source which can feed magnetic fluctuations 
when it exceeds certain thr eshold condit ions. The thresholds can be expressed as A = l—2//3\\ 
for the firehose instability (jParkerlll958l ). the lo wer bound of the temperature anisotropy, an d 



A = 1 + for the mirror mode instability (IChandrasekhar et al.lll958t lHasegawa 



19691) 



the upper bound. The ion cyclotro n instability, com peting with the mirror mode, has the 
onset condition A = 1 + 0.64/ f3^ Al (IGary et al.lll997l ). Here the perpendicular and parallel 
plasma betas are defined as /3± = P±/ (B 2 /2p ) and (3\\ = P\\/(B 2 /2/i ). 



For a perpendicular shock, equations (2) and (3) give the shock strength 
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The temperature anisotropy A 2 can be obtained from equations (4), (5) and (6), 
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(2Ai + l)r 2 r% ' r o 
For 2M\ X > p lt we have (3 ±2 ~ 2M| 1 /r 2 , so the shocked plasma would have a high beta 
independent of j3i, consistent with the findings in the magnetosheaths of outer planets 
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( IRussell et al.l 119901 ). As can be seen from the mirror-mode threshold, high values of (3± 2 
favor the onset of the mirror mode instability; similarly, A 2 ~ l/(r s — 3) under the same 
condition, so A 2 > 1 since the shock strength cannot be larger than 4, also favoring the 
mirror mode onset. 

Figure 1 shows the temperature anisotropy A 2 for r s = 3 and A\ = 1 as a function of 
M/i! and /3i. Only values of < A 2 < 2 are shown. The entropy, S = ^k B \n(P]_P\\/ p 5 ), is 
required to increase across a shock by the second law of thermodynamics; the region which 
breaks the entropy principle is show as "Forbidden" in Figure 1. The majority of the allowed 
area has A 2 > 1 as expected and is preferentially unstable to the mirror mode instability. 
The mirror mode has a lower threshold than the cyclotron mode in this plasma regime, so 
the temperature anisotropy would be quickly reduced by the mirror mode before the ion 
cyclotron mode could develop. The TS with (3\ ~ 32.8 and Max — 16.3 (see §3), indicated 
by the plus sign, is located slightly above the mirror mode threshold. Interplanetary and 
planetary bow shocks have a large variation in Max an d /3i; many of them would also be 
above the mirror mode threshold as indicated by the large area with A 2 >1. Consequently, 
mirror mode instabilities should be a frequent feature downstream of quasi-perpendicular 
shocks. 

For a parallel shock, the magnetic field does not go into the expression of the shock 
strength. The temperature anisotropy 

A _ lfe(^ + |-A) + M li^ + 2 /^-3) nm 
A T + 2M i 1 (l-l/r s ) ' lUj 

and the downstream parallel plasma beta 

/3||2=2^i +2Mll(1_1/rs) - (U) 
For 2M A1 > Pp ~ 2M A1 , which makes the thresholds close to 1, and A 2 ~ (r s - 2)/2 < 
1, giving favorable conditions for the onset of the firehose instability. Figure 2 displays 
the temperature anisotropy A 2 over various upstream conditions for r s = 3 and Ay — 1. 
Compared with Figure 1, higher values of Max at a given (3\ would otherwise lead to smaller 
A 2 . Most of the area shown has A 2 < 1 as expected for a quasi-parallel shock and is unstable 
to the firehose instability. The TS would induce firehose instabilities in the downstream 
plasma if it were quasi-parallel, as indicated by its location in the plot. Many quasi-parallel 
interplanetary and planetary bow shocks will also give rise to firehose instabilities as implied 
by the large area with A 2 < 1. 

Observations show that a quasi-parallel bo w shock becomes unsteady as the Alfven 



Mach number Max exceeds ~ 3 for (3\ ~ 0.5 (IGreenstadt fc Fredrickd 119790 and is of- 



ten associated with large transverse magnetic fluctuations ( IBavassano Cattaneo et al.l 12000 
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Czaykowska et al.ll200ll ). A closer look at Figure 2 gives A 2 ~ 0.97 at Mai = 2 for Pi = 0.5, a 
noticeable thermal anisotropy but not large enough to exceed the firehose onset; at Mai — 3, 
the thermal anisotropy rises above the firehose threshold, generating firehose instabilities 
which will disturb the shock structure. The firehose instability arises from the fast MHD 
mode and produces large- amplitude transverse waves. The results from this simple ca lcula- 



tion agree with hybrid simulations (IKan fc Swiftlll983l ; iKrauss-Varban fc Omidi 



19911 ) . 



3. Magnetic Fluctuations in the Heliosheath 



VI crossed the TS on 2004 December 16 (day 351) at 94 AU and is making the first 
measurements in the heliosheath. The magnetic fluctuations in the heliosheath are charac- 
terize d by a series of depre ssions in the field magnitude which have been called magnetic 
holes faurlaga et aI.ll2Q06al fbh. These fluctuations are similar to those observed downstream 
of quasi-perpendicular in terplanetary and plane t ary bow shocks that have been identified as 
mirror mode structures (Kaufmann et al.l Il970l ; iTsurutani et al.l Il992l ; IViolante et al.l Il995 



Bavassano Cattaneo et al 



199 



The TS is expected to be highly oblique with 6b n ~ 86° at VI, so we use equations (8) 
and (9) to determine whether mirror mode instabilities occur in the heliosheath. MHD 
simulations give the average preshock plasm a density n\ ^ 8 x 10 -4 cm -3 , speed v\ ~ 380 
km s -1 , and temperature T\ ~ 5.4 x 10 5 K (jWhang et al.l 120041 ). which yields Mai — 16.3, 
Pi ~ 32.8 combined with the observed field strength B\ ~ 0.03 nT. The sh ock strength r s 
is ~ 3 estimated from the jump in the fi eld magnitude (jBurlaga et al.ll2005l ) and from the 
spectral slope of TS accelerated particles dStone et al.ll2005l). T he typical thermal anisotropy 
of the solar wind is A ~ 0.7 near the Earth (ILiu et al.ll2006bl ). Expansion of the solar wind 
would further decrease the value of A if the magnetic moment fi ~ T±/B is invariant; 
when the thermal anisotropy exceeds the firehose threshold, firehose instabilities will be 
induced and help suppress the growth of the anisotropy. The two competing processes 
will arrive at a balance close to the threshold, i.e., A ~ 1 — 2//3ii, which gives Ai ~ 0.94. 
Introduction of pickup ions in the outer heliosphere would not significantly change this value. 
The newly-born pickup ions gyrate about the interplanetary magnetic field, forming a ring- 
beam distribution; this configuration is unstable to the generation of MHD waves, which 
scatter the ions quickly toward isotropy (ILee &: Iplll987l ; iBogdan et al.lll99ll ). We have also 



shown that A 2 does not depend on A x when 2M\i ^> f3 x . Substituting the values of Mai, 
Pi, r s and A x into equations (8) and (9) gives A 2 — 1.03 and Pj_ 2 — 42.2. The threshold of 
the mirror mode is 1 + 1/P± 2 — 1.02, smaller than the downstream anisotropy. As indicated 
by Figure 1, even an A 2 slightly larger than 1 can meet the mirror mode onset at high Pi, 
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so mirror mode instabilities should be induced in the heliosheath. 



For the upstream temperature, we have used t he density-weighted average of the pickup 
ion, solar wind proton and electron temperatures (jWhang et al.l 12004 ) . However, the result 
can be shown to be self-consistent. Equation 9 is reduced to fa ~ — z^f! — l/ r s) i n the TS 
case, from which we have the downstream temperature T 2 ~ 2 x 10 6 K given the average 
field magnitude Bi ~ 0.09 nT and density 712 — 0.002 cm~ 3 from equation 7. Approximating 
the TS as a hydrodynamic shock because of the high plasma beta, 

T 2 _ [2 7 M 2 - ( 7 - 1)][( 7 - 1)M 2 + 2] 
T x (7 + 1) 2 M 2 

we obtain T\ ~ 5 x 10 5 K, consistent with the MHD simulation result. Here 7 = 5/3, and 
M ~ 3, the shock Mach number estimated from n%/ 'rt\ ~ 3. 

Given the absence of plasma measurements across the TS, it is hard to estimate the 
uncertainties brought about by the upstream conditions and the shock parameters. A revisit 
to equations (l)-(6) assuming a 10% error in the upstream conditions and the shock geometry 
and strength gives A2 = 1.03 ± 0.42 and a mirror mode threshold 1.02 ± 0.01. Note that the 
uncertainty of A2 is determined largely by the shock strength, since A2 ~ 1/ (r, — 3) in the 
case of the TS. If we use r s = 2.6t_Q^ inferred from the particle spectra (IStone et al.l 120051 ) 
with other parameters fixed, then the temperature anisotropy A 2 = 1.03 - 4.19; the mirror 
mode would be more likely to occur. Interestingly, the shock strength cannot be smaller 
than 2.2 at the present conditions; otherwise A2 would be negative. It should be emphasized 
that turbulence induced by the mirror mode instability would leave the threshold marginally 
satisfied, so we expect A 2 = 1.02 ± 0.01 in the evolved state. 

The mirror mode instability has maximum growth rate at a wave vector highly oblique 
to the background magnetic field. The minimum variance direction of the measured magnetic 
fields has an angle of about 75 ° with respect to the background field when magnetic holes are 
present (IBurlaga et al.ll2006bl ). Mirror mode waves are non-propagating and appear to be 
static structure s, consistent with the observed magnetic fluctuations with a fairly constant 
field direction (IBurlaga et al.ll2006aU bl). Mirror mode waves are also characterized by anti- 
correlated density and magnetic fluctuations. The density measurements are not available, 
so we cannot verify the mirror mode fro m density fluctuations. The density perturbation Sn 
can be estimated from (jHasegawalll969l ) 



Sn 



r 



SB 



n 



T \\ B 

where SB is the perturbation in the field strength. The fluctuation amplitude SB/ B is ~ 0.4 
- 0.9 in the heliosheath, which results in Sn/n ~ 0.01 - 0.03, too small to be resolved by 
future V2 observations. 
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The non-linear evoluti on of mirror mode instabilitie s would also generate holes in the 
background field strength (IKivelson fc Southwoodl Il996l ). The magnetic holes in the he- 
liosheath are of the similar size (in units of the proton gyroradius) to t hose in planetary 



magnetosheaths that are predominate l y produced by mirror mode waves (IKaufmann et al. 



1 970l ; iBavassano Cattaneo et al.l Il998l ; iBurlaga et al.l l2006bl ) . Magnetic holes in planetary 
magnetosheaths have also been explained as slow mode solitons based on Hall-MHD theory 
( IStasiewiczl 120041 ) . The plasma beta in the heliosheath is ~ 40 as estimated above; in this 
case the ion L armor radius is much larger than the ion inertial length, so the Hall-MHD 
theory breaks down (IPokhotelov et al.ll2005l ). 



Deep in the heliosheath, magnetic field lines may be draped and compressed against 
the heliopause if there is no significant reconnection between the solar wind and interstellar 
fields. Analogous to planetary magnetosheaths, the plasma would be squeezed and flow 
along the draped field lines, leading to a plasma depletion layer close to the heliopause. The 
plasma beta is reduced in this layer, so the mirror mode instability may be inhibited and 
ion cyclotron waves may be gene rated. Damping of these waves would suppress the thermal 
anisotropy and heat the plasma (iLiu et al.ll2006bl ). 



4. Summary 

A simple calculation based on temperature anisotropy instabilities explains a variety 
of observations associated with interplanetary and planetary bow shocks. A shock modifies 
the velocity distribution of particles at its surface, producing instabilities downstream of the 
shock which give rise to different types of waves. The calculation also predicts that mirror 
mode waves form downstream of the TS, which is consistent with the observed magnetic 
fluctuations. The present results provide a substantial basis for shock-induced anisotropies 
which act as an energy source for plasma heating in various space and astrophysical envi- 
ronments. 
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Fig. 1. — Temperature anisotropy downstream of a perpendicular shock with r s = 3 and 
Ai = 1 as a function of f3\ and Mai- The color shading denotes the values of A^. The lower 
region is forbidden for a r s = 3 shock since the entropy does not increase across the shock. 
Also shown are the thresholds for the mirror mode (solid line), ion cyclotron (dotted line) 
and firehose (dashed line) instabilities. Regions above the mirror/ion cyclotron threshold are 
unstable to the mirror /ion cyclotron mode and regions below the firehose onset are unstable 
to the firehose instability. The plus sign marks the TS location. 
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Fig. 2. — Temperature anisotropy downstream of a parallel shock with r s = 3 and Ai = 1 as 
a function of f3\ and Mai. Same format as Figure 1. Regions below the mirror/ion cyclotron 
threshold are unstable to the mirror/ion cyclotron mode and regions above the firehose onset 
are unstable to the firehose instability. 



